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Stellar Encounters: A Stimulus for Disc Fragmentation? 
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ABSTRACT 

An interaction between a star-disc system and another star will perturb the disc, possibly 
resulting in a significant modification of the disc structure and its properties. It is still un- 
clear if such an encounter can trigger fragmentation of the disc to form brown dwarfs or gas 
giant planets. This paper details high resolution Smoothed Particle Hydrodynamics (SPH) 
simulations investigating the influence of stellar encounters on disc dynamics. Star-star en- 
counters (where the primary has a self-gravitating, marginally stable protostellar disc, and the 
secondary has no disc) were simulated with various orbital parameters to investigate the re- 
sulting disc structure and dynamics. This work is the first of its kind to incorporate realistic 
radiative transfer techniques to realistically model the resulting thermodynamics. 

The results suggest that the effect of stellar encounters is to prohibit fragmentation - 
compressive and shock heating stabilises the disc, and the radiative cooling is insufficient 
to trigger gravitational instability. The encounter strips the outer regions of the disc (either 
through tidal tails or by capture of matter to form a disc around the secondary), which triggers 
a readjustment of the primary disc to a steeper surface density profile (and a flatter Toomre Q 
profile). The disc around the secondary plays a role in the potential capture of the secondary to 
form a binary. However, this applies only to orbits that are parabolic - hyperbolic encounters 
do not form a secondary disc, and are not captured. 
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1 INTRODUCTION 

Disc fragmentation has been touted as a possi ble process 
for forming brown dw arfs and gas giant planets ( lBosj|l997l : 
IStamatellos et aflboOTh . If the conditions are right, it does indeed 
seem possible that bound objects can be formed in this top-down 
process: what is less clear is if these conditions exist in nature. 
There are two criteria that must be satisfied befor e fragmentation 
can occur, both of which relate to cool ing processes jGammi d200ll : 
iRice et al.ir2003l : iDurisen et alJIlOOTi) . The first condition that the 
disc must satisfy is that it must be (at least locally) gravitation- 
ally unstable under perturbation, to allow fragmentation to begin. 
A disc is gravitationally unstable t o axisymmetri c perturbations if 
the value of the Toomre parameter ( lToomre|[l964) 



Q < 1.5-1.7 for instabilitv (see lDurisen et alj2007h . If this condi- 
tion is satisfied, perturbations can then grow on the local dynamical 
timescale. 

The second condition refers to the growth of these perturba- 
tions. As a perturbation grows, it will undergo heating due to com- 
pression and shocks. Therefore, if any perturbation is to continue 
to grow, it must be able to radiate this energy away efficiently, such 
that the cooling time is comparable to the local dynamic timescale: 
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where is the sound speed, k is the epicyclic frequency (which is 
equal to the angular frequency Q, if the disc is Keplerian), and S 
is the disc surface density. Nonaxisymmetric perturbations require 
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jGammidlloOll : iRice et ai]|2003l : iMeifa et alj|2005t) . If this condi- 
tion is not met, the fragment carmot collapse to a sufficiently bound 
state to survive in the disc (typically, these unbound fragments are 
destroyed by Keplerian shear). Again, cooling processes play an 
important role: if cooling is not efficient enough, the gas will be 
heated by compression and shocks, Cs will increase, and Q will 
therefore increase to above the critical value, stabilising the region. 
Further, most discs will evolve until they reach their minimum Q 
value (representing an equilibrium between heating and cooling in 
the disc). Once a disc reaches this state, it i s said to be marginally 
gravitationally stable i lLodato & Ricel2004h . 
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Recent work (iMatzner & LevinI l2005l: iRafikovl l2005l : 



IWhitworth & StamatellosI l2006l : iRafikovl l2007h suggests that 
the coohng processes in protoplanetary discs are of hmited 
efficacy (which increases the minimum Q value). Recent semi- 
analytical studies of an gular momentum transport in discs by 
lRice& Armitag3 j2009l) and IClarkd ( l2009h show that the inner 
regions of discs are stable against fragmentation, and that gaseous 
material must orbit at distances of at least 20 au (and possibly as 
distant as 70 au) if it is to become susceptible to fragmentation. 
This c onfirms the work of numerical simu lations jStamatellos et al] 
l2007h and is consistent with observations dOreaves et al.l2008h that 
indicate disc fragmentation does not occur in the optically thick 
inner regions of discs, and will only fragment in the outer regions 
of discs with particular properties, e.g. discs that are massive and 
extended. However, there is another way to decrease Q: increase 
E. All the variables in equation (T} (with the obvious exception of 
G) are local: hence it is possible to make the disc locally unstable 
by increasing the local density. This is in principle achievable by 
the presence of another star moving through (or near) the disc. The 
presence of a companion in the star-disc system can draw out the 
disc matter into tidal tails, compressing the gas while reducing 
its optical depth, improving its cooling efficiency. Also, angular 
momentum transfer between the secondary and the disc can trigger 
outward angular momentum transport - this results in inward mass 
transport, increasing density in the inner regions. If this forced 
increase in E by an encounter (which can help satisfy condition 

1) can be accompanied by sufficient cooling (to satisfy condition 

2) , then in principle it is possible that Q can be maintained at 
a low value, and the disc can successfully fragment. Again, this 
illustrates the importance of correctly treating radiative cooling 
in order to correctly ascertain whether disc fragmentation by this 
effect is feasible. 

Some previous simulations suggest encounters can desta- 
bilise disc s, and can cool suffic iently to cause fragmentation (e.g. 
iBoffin et al. 1998; Watkins et al . 1998a,b; Lin et al. 1998), whereas 
others suggest encounters stabilise discs by compressive an d shock 
heating, and prohibit fragmentation (e.g. lLodato et al .l2007h . There 
have been many attempts to simulate such encounters: until re- 
cently, most analyses of the problem have not considered the full 
effects of radiati ve transfer, or hav e used simple cooling time 
parametrisations jLodato et alfcoOTi) . The correct treatment of ra- 
diative transfer (to correctly evaluate the local cooling time) is cru- 
cial: if the disc is compressed by a companion, this compression 
should occur on the dynamical timescale, which is shorter than the 
thermal timescale. Therefore, it is to be expected that the disc's sub- 
sequent evolution will be adiabatic, and hence fragmentation will 
not occur unless the cooling time can be reduced significantly by 
the compression, i.e. 
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This may be possible if the subsequent temperature increase pushes 
the gas into the so-called "opacity gap" (where ice and dust evapo- 
rate), but the decrease in opacity reduces in magnitude as the den- 
sity of t he gas is increased, so fra gmentation still appears to be 
unlikely djohnson & Gammig|2003l) . 

This paper presents the results of a series of high resolution 
Smoothed Pa rticle Hydrodynami cs (SPH) simulations with radia- 
tive transfer jporgan et al] |2009|) of the encounter of a star-disc 
system with a discless companion. The parameter space explored 
by these simulation s has been selected to emulate the work of 
iLodato et all j2007l) , as a means of comparing the differences be- 



tween using a cooling time formalism and more realistic approx- 
imations to radiative transfer. By doing so, the effect of changing 
various orbital properties can be investigated. In summary: no sim- 
ulations carried out in this work exhibit fragmentation. Although 
the encounter does increase E, radiative cooling is insufficient to 
overcome the significant compressive and shock heating. The paper 
is organised as follows: section[2]describes the configuration of the 
numerical simulations performed; section[3]describes the results of 
the simulations, comparing the effects of changing the various or- 
bital parameters; in sections |4]& [5] the results are summarised, and 
the paper is concluded. 



2 METHOD 

2.1 SPH and the Hybrid Radiative Transfer Approximation 

Smoothed Particle Hydr o dynamics (SPH ) fcucvl 1 19771 ; 
lOingold & Monaghad 1 19771 ; iMonaghad 1 19921) is a Lagrangian 
formalism that represents a fluid by a distribution of particles. Each 
particle is assigned a mass, position, internal energy and velocity: 
state variables such as density and p ressure can t hen be calculated 
by interpolation - see reviews by iMonaghanI {iM llool). In 
these simulations, the gas is modelled by 500,000 SPH particles: 
the primary star (and the secondary companion) are represented 
by point mass particles, whic h can accrete the gas particles if 
sufficiently close to be bound dBate et al.lll995l) . 

The S PH code used in t his work is based on the SPH code de- 
veloped bv lBate et af] ( Il995l ). It uses individual particle timesteps, 
and individually variable smoothing lengths hi such that the num- 
ber of nearest neighbours for each particle is 50 ± 20. In the discs 
studied in this work, this provides a minimum disc scale height {H) 
to smoothing length ratio of f» 1 at the disc's inner edge (which 
increases to a maximum of approximately 6 at larger distances from 
the parent star). The code uses a hy brid method of approximate 
radiative transfer jporgan et al]|2009h . which is built on two pre- 
existing radiative algorithms: the polytropic cooling approxima- 
tion d evised by|s tamatellos et al. ( 20 03), and flux -lim ited diffusion 
(e.g. Whitehous e & Batell2004l ; IMaver et alj2007l. see lForgan et al.l 
[2009 for details). This union allows the effects of both global cool- 
ing and radiative transport to be modelled without extra boundary 
conditions. 

The opacity and temperature of the gas is calculated using 
a non-trivial equation of state - this accounts for the effects of 
of H2 dissociation, H" ionisation, He'^ and He+ ionisation, ice 
evaporation, dust sublimation, molecular absor ption, bound-free 
and free-fre e transitio ns and electron scat tering teell & Liiilll994l ; 
iBolev et aL 2007; Sta matellos et al.ll2007l) . Heating of the disc can 
also be achieved by P dV work and by shocks. 
For convenience, a brief recap of one of the test cases undertaken 
by the code is reproduced here . This test case is the evolution of 
a protoplanetary disc, with parameters defined i n a series o f paper s 
by Pickett etal. (2003); Meji'aetal. (2005); B olev et all ( |2006|) ; 
ICaietalJ l l2008l) . The model is a 0.07 Mq Keplerian disc which ex- 
tends to 40 AU, orbiting a star of 0.5 Mq. Initially, the surface den- 
sity profile is E ~ r"^'^^, with a temperature profile of T ~ r^^ . 
The disc is modelled using 2.5 x 10^ SPH particles, with one sink 
particle representing the star. The disc is immersed in a radiation 
field of ro(r) = 3K ; the effects of disc irradiation by the central 
star are not included. 

A snapshot of the disc can be seen in Figure [T] The disc 
matches qualitatively the work of Boley et al, reproducing the 
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Table 1. Summary of the orbital parameters investigated in this work. 



Simulation Mdisc/^Q S oc r ^ M2/M0 Calculated ijperi (au) Actual Rperi (au) e Prograde/Retrograde Inclination 



1 0.1 1 0.1 40 28 1 Pro 0° 

2 0.1 1 0.1 30 25 1 Pro 0° 

3 0.1 1 0.1 50 34 1 Pro 0° 

4 0.1 1 0.1 100 100 1 Pro 0° 

5 0.2 1 0.1 50 36 1 Pro 0° 

6 0.1 1 0.1 50 40 1 Retro 0° 

7 0.1 1 0.1 30 30 1 - 90° 

8 0.1 1 0.1 30 33 7 Pro 0° 

9 0.1 1 0.5 40 10 1 Pro 0° 

10 0.1 1.5 0.1 40 35 1 Pro 0° 




-50 50 

x[AU] 



Figure 1. Snapshot of the lBolev et alj j2006h test disc after 1906 years. 



burst and asymptotic phases, and producing a Q profile consis- 
te nt with their results . A more detailed discussion can be found 
in lForganet^bOOgh . 



2.2 Initial Disc Conditions 

Originally, i t was intended that the disc conditions used by 
iLodatoetalJ ^20071) would be used to facilitate a direct compari- 
son of techniques. However, these initial conditions are scale-free, 
which poses proble ms when radiative transfer is simulated. The 
iLodato et alj bOOTi) parameters were therefore modified to ensure 
that the disc is marginally gravitationally stable. In essence, this re- 
quired decreasing the star mass, and increasing the radial extent of 
the disc. The disc extends from ri„ = 1 au to Tout = 40 au. The 
disc has mass 0.1 Mq, with a central primary star of mass 0.5 Mq. 
The initial surface density profile was chosen to be E oc r~^, with 
a sound speed profile of Ca oc r~* . Two variants of the disc were 
also run: one with a disc mass of 0.2 Mq, and one with E oc r^'^^"^ . 

The discs were evolved in isolation for several Outer Rotation 
Periods (ORPs). This allows the disc to approach an equilibrium 
state and become marginally unstable, and to develop steady-state 
spiral structures (with the exception of the E oc r~^l'^ disc, see 
Figure[2j. As the discs are either stable or marginally unstable, the 
discs will require external stimuli in order to fragment. 



2.3 The Stimulus: Adding a Companion 

With the disc evolved into a quasi-steady state, a secondary star 
was then added (at a separation of 100 au to the primary, in order to 
prevent any non-linear perturbations in the disc by the secondary's 
sudden appearance). The secondary was added with a variety of or- 
bital parameters, comprising a suite of 10 simulations (see Table 
[T]for details). Of the encounters simulated, eight were coplanar, in 
prograde motion; one was coplanar, with retrograde motion; and 
one was inclined at an angle of 90°. Eight simulations were con- 
ducted using the standard 0.1 Mq disc, one was conducted using 
the 0.2 Mq disc, and one using the E oc r~'^^'^ disc. This provides 
data on how disc perturbations vary with the disc mass, surface 
density profile, the secondary's periastron, eccentricity, inclination, 
and the relative angular momenta of the secondary and the disc. 

2.4 Resolution 

It is of utmost importance that any SPH simulations of disc frag- 
mentation correctly resolve the fragments. Bate & Burkert (1997|) 
show that SPH correctly reproduces fragmentation if the minimum 
resolvable mass, 

M,^,„ = 2Nneiahm, = 2Mtot^^^^P^ (4) 

is less than the Jeans mass of interest (m^ is the mass of a single 
SPH particle, Ntot is the total number of particles, Nneigh is the 
typical number of nearest neighbours for a single SPH particle, and 
Mtot is the total mass). The Jeans mass can be estimated using the 
fact that the most unstable wavelength to gravitational instability is 
the disc's scale height, H. This gives a Jeans Mass Mj = Y,H^ — 
Mtot (^)^, and hence 

Mj (hV Ntot _^,lf Ma,a, y Ntot 

Mmin 2 \ RJ Nneigh 2 \ A/, / Nneigh 

All the simulations use 500,000 particles, and have Nneigh ~ 50. 
Using Mmsc = O.IM0, and M* = 0.5Mq, this gives Mj ^ 
200Mntin- With the Jeans mass of interest being two orders of 
magnitude above the minimum resolved mass, this shows that all 
simulations carried out in this work easily satisfy the necessary res- 
olution conditions. 

3 RESULTS 

The results for each simulation are discussed below. A summary 
of all simulations can be found in Table [T] Qualitatively speaking. 
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Figure 2. Snapshots of the three discs used after several ORPs: the 0.1 Mq disc (top), the 0.2 AIq disc (middle) and the E oc r disc (bottom). Note 
each plot has a specific colour bar. 
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Figure 4. The Simulation 1 disc after the encounter. 

the results show features common to all simulations. Therefore, the 
results for simulation 1 (the reference simulation) will be discussed 
in detail - the other simulations will be more briefly described, fo- 
cusing on their unique and differing features. 

3.1 Simulation 1 - The Reference Simulation 

Images of the reference simulation can be seen in Figure [3] the 
azimuthally-averaged radial profiles at the three instants are shown 
in Figures milllTl &[8] As the secondary passes through the disc, 
it imparts significant energy to the disc, causing a strong tempera- 
ture spike at the location of the secondary, and the scale height of 
the disc to be enlarged. The disc exerts significant tidal forces on 
the secondary as it reaches periastron, and angular momentum is 
transferred to the disc. There are three consequences to this: 

(i) The orbital parameters are altered, such that although the ini- 
tial velocity of the secondary is consistent with rperi = 40 au, the 
actual periastron is lower (this can be seen in the E spike in Figure 

13. 

(ii) The secondary can capture disc matter to form a secondary 
disc (M ~ 0.006A/q, comprising around 5% of the initial disc), 
as well as drawing out tidal tails (see Figure[3]l, which depletes the 
primary disc. In total, around 20% of the initial disc mass is lost. 

(iii) The disc is pushed out of equilibrium, and angular momen- 
tum transport is required to stabilise the disc: this results in a steeper 
surface density profile (and a flatter Q profile). This is a common 
feature to all the simulations carried out (although the strength of 
these effects varies with the orbital parameters used). 

Once the encounter is complete, the disc spiral structure is much 
weaker (Figure |4), and the disc has stabilised at all radii, with a 
higher surface density and scale height, and an unchanged temper- 
ature profile. 

3.2 Simulation 2 - A Low Periastron Encounter 

With a low periastron radius in prograde motion, the disc exerts 
stronger tidal forces, resulting in the capture of the secondary on 
a highly eccentric orbit. The secondary captures a similar disc to 
Simulation 1, but the total disc mass lost is slightly higher (due to 
more massive tidal tails being induced). This increased loss results 
in a larger steepening of the surface density profile, as the primary 
disc must make a larger readjustment to retain stability. 
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Figure 5. Surface density profile of the Simulation 1 disc before the en- 
counter (solid line), at periastron (dashed line), and after the encounter (dot- 
dashed line). 




Figure 6. Scale height of the Simulation 1 disc before the encounter (solid 
line), at periastron (dashed line), and after the encounter (dot-dashed line). 
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Figure 7. Midplane temperature profile of the Simulation 1 disc before the 
encounter (solid line), at periastron (dashed line), and after the encounter 
(dot-dashed line). 
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Figure 8. Toorare Q profile of the Simulation 1 disc before the encounter 
(solid line), at periastron (dashed line), and after the encounter (dot-dashed 
line). 

3.3 Simulation 3 - A High Periastron Encounter 

A high periastron encounter, although having the same qualitative 
features as Simulations 1 & 2, has a reduced quantitative effect. 
This is of course consistent with the gravitational influence of an 
object decreasing with separation. The surface density profile is 
steepened, but the effect is only noticeable at large radii, which 
gives a Q profile that remains flat out to around 35 AU (slightly 
further than Simulation 1). 

3.4 Simulation 4 - A Distant Periastron Encounter 

When the periastron is very large, the effects of the encounter are 
reduced significantly. However, this simulation is worth discussing. 
The distance of the secondary from the disc ensures that the disc 
heating is minimal, which allows the outer regions of the disc to 
remain marginally unstable. Also, the surface density profile is 
slightly modified in the iimer regions (Figure |9}, without affect- 
ing the scale height or the temperature profile. This results in a disc 
that is marginally unstable over a larger range of radii (Figure fTOt. 
Although the disc does not fragment in this simulation, it has been 
modified by the encounter to become more amenable to fragmen- 
tation. This suggests that there is a range of periastra for which 
encounters are most effective at destabilising the disc, which will 
depend on the disc's (and the secondary's) properties. 

3.5 Simulation 5 - A Higher Disc Mass Encounter 

The previous simulations have shown that the effect of an encounter 
is to steepen the disc's surface density profile, and to flatten the Q 
profile. Will this behaviour hold for higher disc masses? To inves- 
tigate this, a prograde coplanar encounter was simulated using the 
0.2 Mq disc described in section [Z2l Figure [TT] shows that the 
surface density profile does steepen, but not a great deal: this is 
primarily due to significant mass depletion (and relatively reduced 
mass redistribution) - the secondary captures a disc of 0.01 Mq, 
removing matter from the outer regions to steepen the profile. At 
periastron, a significant increase of the surface density (and scale 
height) can be seen; the increased density at periastron (in com- 
parison with Simulation 1) prevents the temperature increase from 
reaching the same magnitude. The Q profile for the disc (Figure 
I12t was initially rather flat; during the encounter, the companion 
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Figure 9. Surface density profile of the Simulation 4 disc before the en- 
counter (solid line), at periastron (dashed line), and after the encounter (dot- 
dashed line). 
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Figure 10. Toomre Q profile of the Simulation 4 disc before the encounter 
(solid line), at periastron (dashed line), and after the encounter (dot-dashed 
line). 

induces a small region to become very unstable around periastron, 
but due to the mass stripping, the majority of the disc stabilises: the 
region of instability decreases from around 10 — 40 au to around 
10 - 20 au. 



3.6 Simulation 6 - A Retrograde Encounter 

In prograde encounters, the ability of the disc spiral structure to 
couple with the perturber is much improved in comparison with 
the retrograde encounter, and hence angular momentum transfer 
between the perturber and the disc in a retrograde encounter is rel- 
atively smaller (Hall et al. 1996). In general, this is reflected in the 
results obtained. However, depending on the relative phases of the 
spiral density wave and the secondary, the secondary can encour- 
age the spirals to wind more tightly, allowing compressive heating 
in the inner regions. This is evident in the temperature spike at ~ 
10 au in Figure[T3] Despite this, the inner regions of the disc appear 
to be less stable than they were initally. This dependence on the rel- 
ative orbital phases shows that the influence of encounters on disc 
dynamics are indeed more complex than a simple study of orbital 
parameters (such as this work) may initially suggest. 
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Figure 11. Surface density profile of the Simulation 5 disc before the en- 
counter (solid line), at periastron (dashed line), and after the encounter (dot- 
dashed line). 



Figure 14. Toomre Q profile of the Simulation 8 disc before the encounter 
(solid line), at periastron (dashed line), and after the encounter (dot-dashed 
line). 
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Figure 12. Toomre Q profile of the Simulation 5 disc before the encounter 
(solid line), at periastron (dashed line), and after the encounter (dot-dashed 
fine). 
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Figure 13. Midplane temperature profile of the Simulation 6 disc before the 
encounter (solid line), at periastron (dashed line), and after the encounter 
(dot-dashed line). 



3.7 Simulation 7 - A High Inclination Encounter 

As a secondary on a high inclination orbit has a short disc interac- 
tion timescale, it should be expected that the strength of the pertur- 
bation is reduced in comparison to a coplanar orbit. Our results do 
seem to confirm this. Although the outer regions of the disc appear 
to increase in mass (and become less stable), the inner regions of 
the disc are almost unaffected: the Q and S profiles remain similar 
in the region 10 — 30 au, while the temperature profile is unaf- 
fected at virtually all radii. It could be argued that the transport of 
mass outward out of optically thick regions (allowing more efficient 
cooling) is favourable to fragmentation, however none was seen. 



3.8 Simulation 8 - A Hyperbolic Encounter 

As with the high inclination encounter, a hyperbolic encounter has 
a shorter interaction timescale than a parabolic encounter, and so 
should exert less influence on the disc dynamics. Figure [T4l does 
indeed show that the perturbation to Q is weak in comparison with 
Simulation 1: however, the impact of the high- velocity encounter 
has resulted in significant mass stripping, causing a decreased Q in 
the outer regions: however, the velocity dispersion of the outflow 
ensures that the gas is still unbound. 



3.9 Simulation 9 - A High Mass Secondary Encounter 

For the last of the simulations that emulate ILodato et^ ( l2007l) . the 
effect of increasing the secondary mass is studied here. With the 
primary and secondary now of equal mass, the disc begins to feel 
the influence of the encounter much earlier than in previous simu- 
lations. A substantial tidal tail is formed at an early time, and be- 
comes significantly concentrated. This provides significant torques 
at much earlier times than the other simulations, which results in 
a far smaller periastron 10 au), and at greater speed, which in 
turn heats the inner disc significantly, (see Figure |15t . Very effi- 
cient mass-stripping causes the disc to decrease in mass by more 
than half: a large fraction has been swept into the tidal tail at larger 
radii, and into a secondary disc. The significant mass loss and heat- 
ing ensures the primary disc itself cannot fragment - it is now only 
a few au in radius, and extremely hot. 

This simulation would suggest that high mass encounters are 
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Figure 15. Midplane temperature profile of the Simulation 9 disc before the 
encounter (solid line), at periastron (dashed line), and after the encounter 
(dot-dashed line). 
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Figure 16. Surface density profile of the Simulation 10 disc before the en- 
counter (solid line), at periastron (dashed line), and after the encounter (dot- 
dashed line). 



even less likely to promote fragmentation than low mass encoun- 
ters (at least with periastra within the disc). It is possible, however, 
that a distant encounter with a high mass secondary would be more 
successful at producing fragments (cf. Simulation 4), although with 
a high mass secondary's larger perturbative potential, the periastron 
would most likely need to be larger than is considered here, and an 
interesting avenue for further work. 



3.10 Simulation 10 - A Steep Disc Profile Encounter 

Having studied the results of the previous simulations, it seems 
that the result of most encounters is to steepen the surface den- 
sity profile of the disc. This behaviour then begs the question: what 
if the profile is steeper initially? Does the encounter result in an 
even steeper profile, or does there exist an asymptotic profile to 
which a disc will tend (given enough encounters)? To answer this 
question, an additional simulation was run with a disc exhibiting a 
E oc r"''^^ profile. The disc was run in isolation (as with the other 
discs used in this work) for the same timescale: this does mean 
however that the disc is stable over most radii when the encounter 
is begun (although Q is tending towards instability at larger radii). 
Again, this simulation was run for the purpose of studying the ef- 
fects of encounters on the surface density profile, not fragmenta- 
tion, so the need for a marginally stable disc is not so important 
here. 

As can be seen in Figure [16] the effect of the encounter still 
steepens the profile (although comparing with Simulation 1 (Fig- 
ure [S] the magnitude of disc readjustment is indeed smaller). As 
the outer regions of the disc are lower in density than in Simulation 
1, the scale height is more sensitive to the encounter, and there- 
fore see a greater increase as the secondary passes through the disc. 
However, the low density also ensures that compressive heating is 
less effective: therefore a lower temperature increase is seen. 

This additional simulation has shown that surface density 
steepening appears to be common to all discs, regardless of their 
initial profile. However, the steeper the disc initially, the less ef- 
fect the encounter has: this does suggest that there may be a profile 
steep enough that an encounter does not affect it significantly - but, 
this putative profile may be too steep to be physically appropriate. 



4 DISCUSSION 

Having displayed results from a number of simulations, the broad 
trends uncovered by this work will now be discussed. 

4.1 The Influence of Disc Mass, Disc Profile and Secondary 
Mass 

Changing the mass involved in the encounter (and its distribution 
in the system) will of course affect the tidal forces experienced by 
all participants in the encounter. Increasing the disc mass might 
increase its initial instability, but this does not in general indicate 
a more unstable disc after an encounter. The inner regions of a 
more massive disc can remain more gravitationally unstable, but 
the inefficiency of cooling in these regions precludes the possibility 
of fra gmentation ( Ga imniellioOll : iRice et ai]|2003l ; lRafikovll200a: 
iRice & Armitage 200^^ 

Encounters in general steepen the surface density profile: this 
appears to be true regardless of the initial profile, but the magnitude 
of profile steepening does seem to decrease if the initial profile is 
steep. They also flatten the Q profile as a secondary effect, but this 
does not appear to be true if the Q profile is initially flat. 

As far as secondary mass is concerned, the results shown here 
show that its primary influence is increased tidal forces subjected 
on the disc (as well as increased compressive heating at periastron). 
Therefore, for disc-penetrating encounters, a larger secondary has 
a deleterious effect on potential fragmentation in comparison to a 
smaller secondary. 

4.2 The Influence of Periastron Radius 

A crude expectation may be that reducing periastron radius will 
destabilise the disc more, and induce fragmentation more readily. 
This is incorrect, primarily due to radiative effects. The ability of 
the disc to cool has been shown to be an important factor. If the pe- 
riastron is too low, the disc is heated too strongly by the secondary's 
motion through it, and the disc is readjusted into a much more sta- 
ble state due to mass stripping and angular momentum redistribu- 
tion. If the periastron is too high, the disc will not feel the effects 
of the secondary's motion, and the disc will feel no stimulus at all. 
However, if the periastron is within some range of values (deter- 
mined by the properties of the primary, its disc and the secondary). 
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then the result of the encounter is to increase the region of insta- 
bility in the disc. An interesting question that can now be asked 
is: would repeated encounters lead to fragmentation? The answer 
is beyond the scope of this paper, but the results would indicate 
that repeated encounters would in general increase outward angu- 
lar momentum transport (and inward mass transport), pushing the 
surface density profile to an equilibrium value that is too steep for 
gravitational instability to act, as the available mass at large radii is 
insufficient CClark e 2009 ; , Rice & Armitage 200^) . This hypothesis 
requires testing, and is an interesting avenue for further work. 

4.3 The Influence of Angular Momentum Alignment (and 
Inclination) 

As ha s been previously discovered l lHall et alll 19961 : iLodato et"al] 
|2007|) . the most effective encounters are prograde in nature, and 
coplanar to the disc. However, the results of this work show that 
retrograde encounters may be more complex than initially thought. 
The effect of spiral winding by the secondary can encourage com- 
pressive heating in the inner regions, and the angular momentum 
transport in this case allows a disc that is more unstable in the inner 
regions. However, there is no indication that fragmentation would 
be better achieved using anything other than prograde encounters. 

Encounters outside of the disc plane (Simulation 7) can drag 
matter to larger vertical distance, increasing the disc scale height: 
but, the interaction time is insufficient for the encounter to encour- 
age significant change. 

4.4 The Possibility of Binary Capture 

It is difficult to make an encounter precisely parabolic for numerical 
reasons. With the exception of Simulation 8, the encounters simu- 
lated in this work have all been almost parabolic: that is, the total 
energy of the secondary is close to zero, and that any decrease in 
this total energy will result in the secondary becoming more bound 
to the primary-disc system. The efficiency of binary capture can be 
estimated from the change in total energy experienced as a result of 
the encounter, as well as the modification of the orbital eccentric- 
itjQ. Table |2] shows the initial total energy of the secondary (£;), 
the total energy of the secondary after the encounter (Ef), and the 
change in orbital eccentricity (Ae). The total energy decreases by 
a factor of 2 in most cases. This appears to be due to the angular 
momentum transfer between the secondary and the disc it accretes. 
As the secondary leaves the primary-disc system, in order to extract 
material orbiting in the same locale (i.e. to allow it to achieve es- 
cape velocity) it must impart dynamical energy to the gas. Also, the 
process of exciting a tidal tail (from matter at the opposite orbital 
phase) requires angular momentum from the secondary, although 
this appears to be less important. Encounters that do not create a 
sufficient disc (e.g. distant or high-velocity encounters) do not ex- 
perience a significant orbital modification. 

These simulations show that capture efficiency increases with 
decreasing distance, provided that the velocity of the orbit is low 
enough that the secondary captures a non-negligible disc of its own. 
If the orbit is initially strongly unbound (such as Simulation 8), then 
no secondary disc can be captured, and the binary formation effi- 
ciency is effectively zero. Furthermore, a high velocity encounter 

^ These eccentricities are calculated at the initial and final time of each 
simulation assuming two bodies in the system: i) the primary and its disc, 
and ii) the secondary and its disc (if any) - the tidal tail is not considered. 



causes severe disc stripping, which diminishes the possibility of a 
future encounter r esulting in captur e . Thes e results are consistent 
with the work of IClarke & Pringld ( Il99ll) . who study the upper 
limit of capture efficiency for hyperbolic orbits, and find that the 
capture rate in those cases is indeed low. 



5 CONCLUSIONS 

This paper has detailed the results of Smoothed Particle Hy- 
drodynamics (SPH) simulations with hybrid radiative transfer 
jForgan et al1l2009l) of the encounters between a protostellar disc 
system and a secondary star. The influence of different disc pa- 
rameters on the possibility of disc fragmentation was explored, and 
the data is summarised above. In general: e ncounters do not in - 
duce fragmentation (confirming the results of iLodato et al.ll2007l) . 
although there appears to be a subset of orbital parameters that 
modify the disc to make it more unstable over a larger range of 
radii. As to whether this is an indication that there are specific or- 
bital parameters for an encounter that cause fragmentation (e.g. a 
highly eccentric elliptical coplanar orbit which has a distant non- 
penetrating periastron), more work needs to be done. Equally, the 
results presented here show that factors not incorporated in this pa- 
rameter study (such as the relative orbital phases of the primary, 
the disc spiral structure and the secondary) have a non-negligible 
role to play in the resulting disc dynamics, indicating more work 
is required to fully understand their effects. It has been shown that 
angular momentum transfer between the secondary and the disc 
is significant: however, calculations indicate that only around 1 in 
10,000 discs will exper ience encounters of this nature in 1 Myr 
l lClarke&Pringlelll99lh , so this cannot be a typical trigger for an- 
gular momentum transport in discs. 

The most import ant conclusion to dr aw from this work is the 
same as was found bv Lod ato et al.1 ( l2007h - the key parameter that 
determines disc fragmentati on is the cooling rate jGa mmie 200l|; 
iRice et al.ll2003l : [Meiia et ai]|2005l: \slce et alJbOOaT This is inde- 
pendent of whether discs become unstable in isolation, or are driven 
to it by external stimuli such as stellar encounters. 
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